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Abstract 

We consider a model of interacting cosmological const ant /quintessence, where 
dark matter and dark energy behave as, respectively, two coexisting phases of a 
fluid, a thermally excited Bose component and a condensate, respectively. In a 
simple phenomenological model for the dark components interaction we find that 
their energy density evolution is strongly coupled during the universe evolution. 
This feature provides a possible way out for the coincidence problem affecting many 
quintessence models. 



PACS 98.80.-k, 98.80.Cq 



The extensive astrophysical observations of Supernovae (SN) la provide a powerful 
tool to single out the main cosmological energy density contributions at redshift z ~ 1. 
The evidence for a presently accelerated cosmic expansion [1, 2] suggests the presence 
of a new dominant energy source whose equation of state p = up is compatible with a 
cosmological constant. The analysis of the recent experimental data favors in fact the 
range — 1 < u < —0.4 [1, 3, 4]. In particular for a pure cosmological constant the 
favored energy density, in unit of the critical density, would be of the order of f2 A ~ 0.7. 
Unfortunately, this value would be several orders of magnitude smaller than what can 
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be predicted by our present knowledge of fundamental interactions. Moreover, it is also 
quite unnatural that Q\ is just of the same order of the cold dark matter term Q m ~ 0.3 
today, while it has been negligible for, say, red-shifts larger than 1. This is a severe fine 
tuning problem, usually referred to as the coincidence problem. 

A much more appealing scenario is to look at this additional term as due to a scalar 
field whose dynamics makes it negligible at high redshift and dominating the cosmic 
energy density at z < 1 [5, 6, 7, 8]. However this wide class of model is not free from fine 
tuning problem, due to the arbitrariness in the choice of the scalar field potential, but it 
has nevertheless more appealing features. This is for example the case of tracker fields 
scenarios, where the scalar dynamics shows a pseudo-attractor solution, valid for a wide 
range of initial conditions and providing values of and [5, 6] in good agreement 
with SN la observations. 

Recently interacting quintessence scenarios have been extensively studied, considering 
a scalar field coupled with gravity, ordinary matter and dark matter [9, 10, 11, 12, 13, 14]. 
In particular, the recent data on CMBR anisotropy have provided a chance to constrain 
quintessence models [15, 16, 17, 18, 19, 20, 21, 22], though the high degree of arbitrariness 
on the scalar dynamics makes these constraints very model dependent. In this letter we 
assume a more phenomenological scenario, in which the quintessence and dark matter 
components are looked upon as two different components of the same fluid, a condensate 
contribution and a thermally Bose distributed "gas", with energy densities pa and p m , 
respectively. We start considering the case where the condensate energy term corresponds 
to a pure cosmological constant, tv — — 1. As it will be clear in the following, the main 
features of the model still hold for a more general equation of state with uo < 0. The gas 
component p m is assumed to be cold matter, with vanishing pressure. Finally we consider 
for simplicity a spatially flat universe, Qk = 0. 

We parametrize the interaction in terms of powers of the two energy densities p\ and 
p m . The time evolution of radiation energy density, pn, consensate and dark matter are 
then dictated by the following equations 

Pr = -4-p r , (1) 

a 

p m = -3-p m -kM 5 - 4a -^p a m pi , (2) 

(X 

p A = kM*-^p« m p{ , (3) 

where k is a dimensionless coupling constant and M a suitable mass parameter. In terms 
of the redshift z equations (l)-(3) can be conveniently rewritten as 

d -¥ - ih^ ■ « 
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dz 1 + z V y/pr(z)+Pm(z)+p A (z)* 



1 + 2 \/pr(^) + Pm(^) + P~a(z) 



(6) 



where p"j = Pi/p co , with p Co = ^H^MpJ [Sti) = h 2 0.81-10 46 GeV 4 the present critical 
energy density, and 

_ , , f^r"- 3/2 • (7) 

7T V M ) VMV 1 ; 

By defining x = log 10 (1 + z) we then get 

dp r (x) 




dx 



= 4 1og(10)p r (x) , 



d ~ Pm{x) = log(10) L~ Pm {x) + \ , PJM£M I , (9) 



da; 



^pa(^) = _ A log(10) af^f) . (io) 



dx \/ Pr(x) + p m (x) + p A (x) 



Note that for the parameter A < 0, the interaction term on the r.h.s. of (9) and (10) de- 
scribes a continuous transfer of energy from the condensate to the ordinary non-relativistic 
particles. The corresponding scenario is a universe filled by both the condensate and a 
normal gas of dark particles, but with an effective mass of the condensate larger than the 
mass of thermal excitation in the gas phase. Direct decay of the condensate in dark mat- 
ter particles is expected to be proportional to p A only, while stimulated decay mechanism 
is of course proportional to number density of dark matter particles too. The natural 
expectation in this case is an interaction term with a — /3 — 1. In the following we will 
keep varying these parameters, to study the general features of the model as function of 
a and j3 and eventually study in more details the case whereboth are set equal to unity. 

We will be mainly interested to the case A < 0. The case A > 0, which will be 
briefly treated in the following, also represents a possible framework. In this case, as 
the temperature lowers, a fraction of the energy stored in the gas of dark particles flows 
in the condensate, as during a temperature driven phase transition. Unfortunately this 
model suffers of the very same coincidence problem of many other quintessence scenario 
and thus it is much less appealing. 

The system (8)-(10) can be numerically solved, as a function of x, by fixing the values 
of p r , p m and p A at present (#j = 0), and then following the values of the energy densities 
back in time. In particular, to describe the features of this coupled system, we have chosen 
as reference values p m (0) = Q m = 0.3, pa(0) = Qa = 0.7 and the equivalence red-shift 
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z eq ~ 3100. This implies 

p r (z) p r (0) p r (0) 1 n , 

Pm(*) P m (0) 1 + J An(0) 1 + ^ ' 1 ^ 

It may seem that fixing the values of the energy densities parameters today is already 
a fine-tuning condition imposed to the evolution. This of course would be true if no 
interactions would be considered in the model. We will show in fact that, since p\ and 
p m evolution are coupled, the fact that their ratio Q\/Q m is of order one today represents 
a much less ad hoc condition for a wide range of the parameters a, (3 and A. The fact 
that Q\ ~ Q m represent the main contribution to the energy density today is ultimately 
related to the fact that the interaction term keeps their evolution quite close throughout 
the universe expansion. We can say that the only remaining fine tuned condition is given 
by the value of the equivalence redshift, where radiation start becoming negligible. 

In order to study the role of the free parameter a, 7 on the evolution of the energy 
densities, we start fixing for example A = — 10~ 3 and choose a = 1 and (3 = 0. In Figure 
1 are reported the energy densities for radiation p r , matter p m and condensate p\, as a 
function of x. For these values of free parameters, the interaction is due to p m only, but 
differently from the noninteraction case, the energy density of condensate here evolves 
in time. The evolution is of course enhanced for larger values of a, but for x < 2 the 
condensate energy density p\ reaches a constant value which, at present, provides an 
effective cosmological constant. Thus for (3 = the model predicts scenarios very similar 
to other quintessence models and no new features overcoming the coincidence problem 
appear. 

The behavior of the system completely changes if one takes non vanishing values for 
(3. In this case, by increasing the value of (3, the values of dark energy p\ and dark matter 
p m are more and more connected. As it is clear from Figure 2, where the energy densities 
are plotted versus x for different values of (3 and for a = 1 and A = —0.001, due to the 
interaction term Eq.s (9) and (10), p m decreases for large redshift with the increasing of 
(3, whereas p\ still remains with larger values. As a consequence their evolution remain 
closer. The two quantities start to intersect each other for values of (3 larger than 0.8. 
This is a relevant new feature of the model which solves the problem of coincidence, 
simply because the situation where the two densities p A and p m are of the same order of 
magnitude it is not a peculiarity of the present epoch but occurred many times in past. 

In Figure 3 we show the plots corresponding to the non-interacting case A = and the 
evolution for A = —0.001 and for a, (3 = 1 . In this case for < x < 2 the two graphs 
exactly overlap, since in Eq. (9) the first term dominates whereas in Eq. (10) the l.h.s. 
is almost vanishing. This feature of the model is extremely important because provides 
a proper behavior of dark matter and energy for redshift smaller than z eq which does not 
affect large scale structures formation. 

The behavior shown by the energy densities plots for A = —0.001 and (3 = 1, reported 
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in Figures 2 and 3 can be easily understood. By increasing the redshift x the adimension- 
alized Hubble parameter H = ^/p r + p m + p\ at the denominator of Eq.s (9) and (10) 
also increases, but less than the product of energy densities at the numerator does. Thus 
for large values of p m the interaction term in Eq. (6) starts to be effective, consequently 
Pa increases hence enhancing the interaction term. This occurs till p\ reaches the dark 
matter contribution (at x ~ 3 for this choice of parameters). From this value on p m feels 
the presence of p~A and starts to decrease whereas the dark energy remains constant. This 
behavior ends when the first term in the r.h.s of (9) becomes the dominant one and thus 
p m starts to increase again and so the cycle repeats. 

In Figures 4(a) and 4(b) we report the quantity log 10 (p m /p A ) versus x for (3 = 0.6 and 
P = 0.9, respectively. The other parameters are fixed to be a = 1 and A — 0.001 as for the 
plots of Figure 2. As it is clear from Figures 4 for (3 = 0.6 the interaction term is too small 
to produce sizeable effects. This situation drastically changes for (3 — 1 where the ratio 
oscillates with peaks and dips corresponding to the cycles of the Figure 2. Actually the 
value (3 ~ 1 is the value which most naturally fit the stimulated decay scenario previously 
outlined. 

The peculiar features of the model do not critically depend on the value of A. In Figure 
5 the energy density plots are shown for several values of A in the range — 1 < A < — 10~ 9 . 
As it is clear from the figure, the redshift at the equivalence increases for smaller value of 
|A|. If one imposes the equivalence redshift in the interval z eq = 3100^400 [23] then A must 
satisfy the bound |A| < 10~ 4 . Reminding the definition of A (see Eq.(7)), it is possible 
to recast the bound on this parameter, for a given choice of a and f3, in a bound for the 
mass parameter M and the dimensionless coupling k defined in Eq.(7). In particular for 
a — (3 — 1 one would get 



For values of the mass scale M already larger than, say, few GeV, Eq.(12) is not con- 
straining at all the coupling k. 

The main result of SN la data is that expansion is accelerated for z ~ 1. In order to 
understand if our model is indeed compatible with this result we consider the effective 
equation of state of the primordial fluid w(x), defined as 

,x = Pr(x)+p m (x) +p A (x) = e3&-p A (x) 

~ p r (x) + p m (x) + p A (x) p M {x) + Pr{x) + p A (x) 

In particular a phase of accelerated expansion is guaranteed if w < — 1/3. It is also worth 
noticing that recent data on the supernova SN1997ff at z = 1.7 indicates, moreover, that 
the universe expansion for this redshift is still decelerated. This is actually in beautiful 
agreement with our results, as can be seen in Figure 6, at least for this choice of parame- 
ters. Figures 6, and 8 and 9 also shows that the accelerated phase show that accelaration 
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occurs for z < 1 without assuming any particular fine-tuning of free parameters. The 
behavior of w(x) for larger redshifts and for two different values of A is reported in Figure 
7. The deep occurring for x ~ 3.5 is a genuine effect of the interaction terms, which is 
more evident for large values of f3 and/or of A as shown in Figures 8 and 9, respectively. 
Note that the deep appears for |A| > 10~ 5 , and with the increasing of |A| moves to smaller 
values of x. 

Concerning the prediction for the age of the universe, as function of p~\ and p m , this 
does not provide any severe bound on the free parameter ranges. For all choices of the 
parameters which satisfy z eq = 3100^400 we have checked that the result is the same, in 
the present uncertainty range, of a customary non-interacting cosmological constant. This 
is simply due to the fact that, during the last few red-shifts, which mainly contribute to 
the total age of the universe, p\ and p m effectively behave as a cosmological constant and 
non interacting cold matter. 

The choice u>\ = —1, even if affects of course the time evolution of dark matter and 
dark energy, is not crucial to get the main features of p\ and p m evolution. In particular 
a different choice for u\, provided that uj < —0.4 is still producing the behaviour already 
discussed , with p\ and p m following and intersecting each other during the universe 
evolution. This is shown in Figure 10, where the energy densities versus x are plotted for 
a — P — 1, A = —0.001 and uo\ = —0.5 and —0.8 respectively. 

As stated above the class of models with A > also describe a possibly interesting 
scenario. In this case, in fact, the condensate is not produced till the redshift (temper- 
ature) does not goes under a critical value, x c . This is shown in Figure 11 where for 
the corresponding choice of parameters x c ~ 2.8. For x larger than x c , p^ is vanishing, 
and p m and p r have the same behavior they would have in absence of interaction. At 
x = x c the fluid undergoes a phase transition and part of the energy density stored in 
dark matter flows in the condensate till it start dominate at low red-shifts. Unfortunately 
this scenario is still affected by the coincidence problem. The value of x c where the phase 
transition takes place is determined by the values of p m and p\ today or, viceversa, it 
should be accurately tuned to provide the observed values for the matter /cosmological 
constant energy ratio today. 

To summarize, the coupled quintessence framework provides an interesting scenario 
where the mutual interaction between dark matter and dark energy components shows 
a quite close evolution of their energy densities during the universe expansion. This 
remarkable feature represents a simple way out to the coincidence problem, since it makes 
the present situation of p\ ~ p rn a cyclic condition which have occurred many times in 
the past. The detail of the energy density evolution is of course a function of the assumed 
interaction term. We have considered a simple model where the dark energy is provided 
by a condensate component of a fluid. The interaction with the gas of non relativistic 
particle (dark matter) is then due to stimulated decay processes of the condensate into 
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gas particles. Though the model is still lacking a detailed microscopic description, it is 
nevertheless quite phenomenologically well motivated, and it would be worth analyzing 
their properties in more details, in view of the natural way it deals with the cosmic 
coincidence problem. 
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Figure 1: Energy densities in unit of p co for a = 1, (5 = and A = —0.001. 
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Figure 2: Energy densities in unit of p Co for a = 1 and A = —0.001 and varying /?; the 
corresponding value of (3 is reported on the top of each graph. 
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Figure 4: The log 10 [p m / Ph\ versus x for the two values (3 = 0.6 (Fig. (a)) and (3 = 0.9 
(Fig.(b)). The values of a and A are the same of Figure 2. 
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Figure 5: Energy densities in unit of p co for varying A; the corresponding value of A is 
reported on the top of each graph. The quantities a and (3 are both fixed to be equal to 
1. 
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Figure 6: Total equation of state for A = —0.001 and a = (3 = 1. The behaviour of w is 
plotted for a range of z near the present epoch. 
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Figure 7: Total equation of state for A = (w ) and A = — 0.001 (w) with a = (3 — 1. 
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Figure 8: Total equation of state for different values of (3, with a = 1 and A = —0.001. 
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Figure 9: Total equation of state for different values of A, for a = (3 = 1. 
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Figure 10: Energy densities in unit of p Co for a = (3 = 1, A = —0.001 and uja = —0.5 and 
—0.8 respectively. 




Figure 11: Energy densities in critical units for A = 0.0014 and a = (3 = 1. 



15 



